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Abstract 

We demonstrate how some problems arising in simplicial quantum gravity can be suc- 
cessfully addressed within the framework of combinatorial group theory. In particular, 
we argue that the number of simplicial 3-manifolds having a fixed homology type grows 
exponentially with the number of tetrahedra they are made of. We propose a model of 3D 
gravity interacting with scalar fermions, some restriction of which gives the 2-dimensional 
self-avoiding-loop-gas matrix model. We propose a qualitative picture of the phase struc- 
ture of 3D simplicial gravity compatible with the numerical experiments and available 
analytical results. 



1 Introduction 



The success of the matrix models as a theory of 2D quantum gravity and non-critical 
strings has brought about a hope that analogous discrete approaches might be instructive 
in higher dimensions as well. The most natural way to introduce discretized quantum 
gravity is to consider simplicial complexes instead of continuous manifolds. Then, a path 
integral over metrics (which is the crux of any approach to quantum gravity) can be 
simply defined as a sum over all complexes having some fixed properties. For example, 
it is natural to restrict their topology. In the present paper, we consider only the 3- 
dimensional case, where, by definition, a simplicial complex is a collection of tetrahedra 
glued along their faces in such a way that any two of them can have at most one triangle 
in common. 

To introduce metric properties, one can assume that complexes represent piece-wise 
linear manifolds glued of equilateral tetrahedra 0. The volume is proportional to the 
number of them. In the continuum limit, as this number grows, the edge length, a, simul- 
taneously tends to zero: a — > 0. In three dimensions, there is a one-to-one correspondence 
between classes of piece-wise linear and smooth manifolds. In other words, any homeo- 
morphism can equally be approximated by either a piece-wise linear or smooth map. This 
property makes the model self-consistent. 

Within a piece-wise linear approximation, the curvature is singular: the space is flat 
everywhere off edges. Therefore, one should consider integrated quantities, e.g., the mean 
scalar curvature 

r 1 
J d 3 x y/gR = a(2vriVi - 6iV 3 arccos -) (1.1) 

where a is the lattice spacing and No, N±, jV 2 , N 3 are the numbers of vertices, edges, 
triangles and tetrahedra, respectively. 

For manifolds in 3 dimensions, the Euler character vanishes \ = Nq — Ni+N 2 — N% = 0. 
Together with the other constraint N 2 = 2 A3, it means that only 2 of the AVs are 
independent, Aq = Nq + A3, and a natural lattice action depends on 2 dimensionless 
parameters. Now, we are in a position to define the simplicial-gravity partition function 

B S §§: 

Z(a,/i) = J2e aN °-^ = J2Z N ^)e~" N3 (1.2) 

{C} N 3 

where S{c} is a sum over some class of complexes (e.g., of a fixed topology). 

It can be shown that, if A/3 is fixed, Nq is restricted from above as No < ^N^ + const. 
Therefore, it is natural, keeping a fixed, to tend fi to its critical value, /i c , at which the 
sum over A/3 in Eq. ( |1 . 2|) becomes divergent. In the vicinity of /i c one could expect to 
find critical behavior corresponding to a continuum limit of the model (of course, it can 
require a fine-tuning of the second parameter a). 
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For the model to be physically consistent, the free energy per unit volume log Z^, A (a) 
has to be finite, which means that the series over N3 in Eq. ( |1.2| ) has to have a finite radius 
of convergence: < /i c < +00. If it holds, statistical models defined on such an ensemble 
of fluctuating lattices possess reasonable thermodynamical behavior. Thus one can hope 
to describe matter coupled to quantum gravity. This program has been successfully carried 
out in two dimensions with help of the matrix model technique [0, §]. In the 3D case, 
analytic tools are lacking and numerical simulations have so far been the only source of 
information [||. Q. 

Many questions, which are simple in two dimensions, become stubborn in three. For 
example, 2-manifolds are classified according to values of the Euler character, therefore 
topology, homotopy and homology are actually equivalent descriptions. In three dimen- 
sions, sets of manifolds characterized according to these 3 criteria are essentially different. 
A natural question to ask is how large a class of complexes {C} in Eq. ( |1.2| ) can be. We 
shall argue that it is sufficient to fixe only the first homology group Hi (C, Z) in order to 
have a finite value of /i c . It can be reformulated as the following statement: "The number 
of simplicial manifolds constructed of the given number of tetrahedra, A3, and having a 
fixed homology type (viz., the first Betti number and the torsion coefficients) grows at 
most exponentially function of A3". 

Let us notice that the existence of the exponential bound is not broken down, if one 
weights every complex with a positive weight growing at most exponentially with A3. As 
well, if one manages to represent a sum over complexes as a weighted sum over another 
class of objects and to prove that the weights are exponentially bounded, then one can 
simply estimate the number of objects in this new class ignoring the weights. 

The simplest example of a weighted sum of complexes is given by the canonical parti- 
tion function Z^ 3 (a) introduced in Eq. ( |1.2| ). Another example of physical interest is the 
partition function in the presence of free matter fields. To introduce it let us consider the 
set of polyhedral complexes dual to simplicial ones. Their 1-skeletons are some 4-valent 
graphs, whose vertices correspond to tetrahedra, and links are dual to triangles. 

Such a graph can be defined by the adjacency matrix 



^(4) I 1 if vertices i and j are connected by a link _ ^. 

v I otherwise 



Of course, the set of 4-valent graphs is not identical to the totality of simplicial com- 
plexes. Similarly, in two dimensions, the set of ordinary <f) 3 Feynmann diagrams is different 
from that of triangulations - one has to introduce the notion of fat graphs to establish 
the equivalence. 

If one attaches the n-component vector xf (// = l,...,n) to the z-th vertex, the 
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corresponding gaussian integral can be performed explicitly for every given graph G^': 



r N 3 -l n -i N 3 -l n 

Z matt er{G^) = J n n^I^E G S } « - = (det'L^p (1.4) 



. JV3-I n ^ N3-I 

II II dxfexp [- - £ 
i=i /^=i *j=i 

where the discrete Laplacian is given by 



L% = 4% - Gf (1.5) 

There is the nice combinatorial representation of the determinant given by Kirchhoff's 
theorem PI: 



det'L = \T(G)\ (1.6) 

where T(G) is the set of all possible spanning trees embedded into a graph G or, equiv- 
alently, a number of connected trees which can be obtained from the graph by cutting 
its links. To kill the zero mode in Eq. ( |1.4j ), we fixed the field at the iVyth vertex. The 
theorem states that det'L does not depend on a choice of the vertex. Therefore, |T((j)| 
is the number of rooted trees. 

The number of spanning trees of a fc-valent graph Gffi with n vertices can be estimated 
from above as Hi 



The case we are particularly interested in is the 2-component Grassmann field, where 
we obtain the gravity+matter partition function in the form 

Z«(a,/i) = ^ det'L e aNo ~^ = £ zg(a) e -^ 3 (1.8) 

{C} N 3 

The weights, det 'L, are positive integers, hence, 

Z N3 (a)<Z§(a) (1.9) 

In two dimensions, this type of matter has the central charge c = —2 and the corre- 
sponding matrix model has been solved explicitly by purely combinatorial means [|7| . We 
can repeat the same trick in 3 dimensions: 



Z«(a,/i) = ^det'Le^ "^ 3 = £ ]T e aN °~^ = £ £ e «"°-^» (1.10) 

{C} {C} {T(C)} {T} {C*(T)} 
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Namely, we have the sum over complexes, J2{c}i weighted with the number of spanning 
trees in the dual 1-skeleton of each complex, Yl,rf(c)}' ^ e can obtain exactly the same 
configurations by taking the sum over all possible trees, XVn.; weighted with the number 
of complexes, L C mp which can be recovered from a given tree T by restoring cut links. 

In three dimensions, the appearing trees correspond to spherical simplicial balls ob- 
tained starting with a single tetrahedron by subsequently gluing other tetrahedra to faces 
of the boundary. If we denote as n , rii and n 2 numbers of vertices, edges and triangles 
on the boundary of a ball, we find that 



n = N 3 + 3 n 2 = 2(N 3 + l) (1.11) 

In contrast to the 2-dimensional case, it is extremely difficult in three dimensions to 
control a topology of complexes constructed from a ball by pairwise identifications of 
triangles belonging to the boundary. Therefore, we shall not address the homeomorphic 
properties of complexes in this paper. On the other hand, homotopy information is 
nicely coded in the fundamental group, which is going to be the basic object for us. The 
combinatorial group theory provides us with the most adequate mathematical framework 
within which all physical questions can be asked and, hopefully, answered. 

The outline of the paper is the following. 
In Section 2 we collect some basic definitions and results from combinatorial group theory 
and 3-dimensional topology. 

In Section 3 we give an interpretation of the 3D gravity partition function in terms of 
group presentations. 

In Section 4 we argue that the number of simplicial complexes of a given volume having 

fixed homology type grows exponentially as a function of the volume. 

In Section 5 we introduce a reduced model and establish a connection with the 2D loop 

gas matrix model. 

Section 6 is devoted to a discussion. 

The present paper can be regarded as an extended version of the previous one 0. A 
reader must be aware that Ref. || contains some incorrect statements. In particular, the 
proof of the main theorem stated in it is not mathematically satisfactory. 



2 Complexes and group presentations 

Although almost all the mathematical material collected in this section can be found in the 
standard reference books on lower dimensional topology and combinatorial group theory 
(see for example [|K], [ll]]), it is convenient to repeat it here for reader's conveninience. Let 
us start with some definitions. 

If it is not specifically stated, by a complex we always mean a finite cell (polyhe- 
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dral) decomposition of a compact closed orientable 3-manifold. Any complex, C, can be 
constructed by first assembling the 1-dimensional cells (giving the 1-skeleton, Ki), then 
attaching 2-cells (giving the 2-skeleton, K 2 ) and analogously for K 3 . 

A simplicial complex is a complex in which /c-dimensional cells are /c-simplexes, i.e., 
points, segments, triangles and tetrahedra. The intersection of any two simplexes either 
is empty or consists of exactly one less dimensional simplex. Any simplex in a simplicial 
complex can be unambiguously determined, including an orientation, by a list of vertices 
from its boundary. It means, in particular, that any edge connects a pair of distinct 
vertices. 

The fundamental group tti(C) is determined by the 2-skeleton of a complex: ni(C) = 
iTi(K 2 ). The standard algorithm to read off a presentation 7r 1 (C) is the following: 

1. Construct a spanning tree T of the 1-skeleton of C. 

2. Put into correspondence a generator, Oj, to every link of K\ which is not in T (i.e., 
to each cut link). Links belonging to T are formally associated with the identity, I. 

3. Fixe some orientation of the links of K\\ a change of the orientation corresponds to 
the inversion of the generator: a, — > a^ 1 . 

4. For every face in the 2-skeleton, the cyclic order of the oriented links forming its 
perimeter gives a relation: rj = a^a^ . . . = 1; q = ±1 according to the orienta- 
tion of the i'th link (one puts the identity, 1, on places corresponding to the links 
which are in T). 

As any spanning tree cuts all loops in Ki, all relations obtained in this way are non-trivial 
(the trivial relation is, by definition, of the form 1=1). 

In this way one finds the fundamental group presentation 



where m equals the number of 2-cells in C and n is not bigger that the number of 1-cells. 

This procedure gives a nonreduced group presentation in the sense that presentations 
differing by the trivial cancellations, acr 1 — > 1 and a~ 1 a — > 1, correspond to distinct 
complexes. The simplest example is the following identical presentations of the trivial 
group 



In both cases K 2 consists of the single triangle but, in the second case, its edges are 
identified in such a way that the corresponding complex has only one 0-cell and only one 



7Ti(C) = (a x a 2 • • • a n \r\r 2 ...r m ) 




P' E = (a | a, a, a) 
P'e = 



(2.2) 
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1-cell. The 2-skeleton here is the so-called Zeeman's "Dunce hat"^, a closed contractible 
fake-surface. Both 3-complexes represent S 3 JT2]. 

Therefore, in order to maintain correspondence between presentations and complexes, 
one has to consider nonreduced presentations. However, as we shall see later, in the case 
of simplicial 3-complexes, this subtlety can easily be overcome. Therefore, in what follows, 
we shall use the term "presentation" hoping it will always be clear from a context what 
is meant in every particular case. 

Two presentations describe the same complex, if they can be transformed into each 
other by a sequence of the following transformations: 

1. permutation of generators or relators; 

2. inversion of a generator, relator, r = a\di . . . — > r 1 = a k x . . . a 2 1 a 1 1 ; 

3. cyclic shift of a relator: r = aia2 . . . a*. — > afVai = a 2 . . . a^ai; 

These moves obviously divide the totality of all presentations into equivalence classes. We 
shall always identify presentations from the same class. 

It is known that all 3-manifolds possess finitely presentable fundamental groups. 
Therefore, any of them can be described as the factor group 

G = F(X)/N (2.3) 

where F(X) is the free group on n generators, X = ■ ■ -CL n ), and N is the normal 
closure in F(X) of a subgroup generated by some finite set of words R C F(X) . The words 
constituting R are the relators and this is an equivalent description of the presentation. 
R form a basis of the subgroup (R). Any other set of generators of (R) can be obtained 
from R by applying the Nielsen transformations: 

Nl) replace rj by rj 1 . 

N2) replace rj by r^i with i ^ j. 

This pair of transformations generate a group including permutations as a subgroup. Any 
sequence of the Nielsen transformations carries a presentation (X\R) into an equivalent 
one. Let us notice that they do not change the number of generators. 

To obtain the normal closure of (R) , one has to add the third move, namely, the con- 
jugation by an arbitrary word w 6 F(X): 



*' I am indebted to J.Howie for pointing me out this example. 
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N3) replace Tj by wrjw \ 

In order to establish the isomorphism of two groups, one has to use the Tietze trans- 
formations: 

Tl) add or remove a relator which is a consequence of other relators: 

(X\R) ^ (X\R,y) 
where y € N, N is the normal closure of (R) in F(X). 
T2) add or remove a generator which is a word in the others: 

(X\R) < — ► (X,x\R,x- y w) 
where R and w do not involve x; w is arbitrary. This move is called the prolongation. 

It can be shown that these moves connect any two presentations of the same group. 
It should be said that there is no correspondence between the Tietze moves and the 
transformations of complexes: in most cases they yield presentations having no realization 
in the form of 3-complexes. 

One of the most profound connections between homotopy theory and group presenta- 
tions is established via the so-called formal defomations of complexes ||13|| . Let us introduce 
this notion. An elementary n-expansion, C /" C", is an attachment to C of an n-ball 
B n along all of the boundary dB n except one (n — l)-cell: C" = C'\J B n , C f]B n = B n _i. 
If the complexes in question are simplicial, the ball is simply the single tetrahedron. An 
elementary n-collapse, C \ C", is the inverse of the elementary n-expansion. A formal 

n 

n-deformation, C /\ C", is a finite sequence of the elementary expansions and collapses 
in which the maximum dimension of the balls equals n. 

If a complex C can be n-deformed to a point, *, one says that C is n-collapsible. This 
term should not be confused with contractability, which stays for the simply-connected 
complexes, vr 1 (C) = 1. 



As was proven by P.Wright [13], any two contractible 2-complexes, Tri(K' 2 ) = iti^K'^) 



1, can be transformed into each other by a finite sequence of the formal deformations not 

3 

exceeding the dimension 3, K' 2 /*\ K 2 \ if and only if the corresponding presentations 
read off in the standard way from the two 2-complexes K' 2 and K 2 can be transformed 
into each other by a finite sequence of the Nielsen transformations extended by the trivial 
prolongation. We shall call this set of moves the AC transformations: 

AC1) cancellation of aa" 1 and its inverse. 
AC2) ti — > r^ 1 and — > r^j with i ^ j. 
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AC3) Ti — > wrjW 1 where u> G -F(-X'). 

AC4) (X\R) < — ► (X, a\R, a) where a£X,a£R. 

These transformations will play an important role in the sequel. 

It is not known whether arbitrary two 2-complexes having the same fundamental group 
can always be connected by a chain of the formal 3-deformations, in other words, whether 
in three dimensions the simple homotopy type and the 3-deformation type coincide. The 

3 

famous Andrews- Curt is conjecture reads: for any contractible compact Ki, K2 /*\ * 



holds. It seems that this conjecture is commonly believed to be false |15[]. In any case, as 
we have decided to forget about the notion of homeomorphy, the 3-deformation classes of 
2-spines of simplicial 3-manifolds seem to be the most natural finest classification in our 
context. 

2 

Another theorem of P.Wright jnj] says that K' 2 /*\ K% holds if and only if the cor- 
responding presentations can be transformed into each other by a sequence of arbitrary 
prolongations (T2). It is important that the cancellation of aa" 1 is not allowed as an 
independent move. For example, Zeeman's "Dunce hat" is not 2-collapsible albeit con- 
tractible. 

It is known that the first homology group is the abelianized fundamental group: 
Hi(C,Z) = 7r 1 (C)/[7r 1 (C), 7i"i(C)]. In the abelian case, it is convenient to represent the 
defining relations in the form of a linear system of equations: 

n 

P= (xi...x n \Y,RijXj) ( 2 - 4 ) 

j=l 

where the relation matrix R has integer entries. 

Another set of free generators can be obtained by rotating the vector (x\ . . . x n ) by an 
invertible unimodular integral matrix: 

n 

x >. = M ijXj M e SL(n, Z) (2.5) 
i=i 

and analogously for relators. 

It is known that pre- and post-multiplications by SL(n, Z) matrices can be expanded 
into sequences of the following elementary operations 

1. permutation of rows. 

2. multiplication of all elements in a row by —1. 

3. addition of a row to another row. 
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and the same for columns. 

These operations correspond to the Nielsen transformations and give a homomor- 
phism of the group of the Nielsen moves into SL(n, Z). The kernel is generated by the 
conjugations, rj — > r^jr^ 1 . 

An arbitrary rectangular n x m matrix can be transformed into the diagonal form 

D = diag(di,d,2, . . . ,dk) k = mm(n,m) (2-6) 
where di are non- negative integers and di divides d i+ i (di\d i+1 ). 

Thus, any abelian group has the unique presentation of the form 

G — (xi . . . | x-^ 1 , x^ , • • • , x ' q , — Zp^ x Zp 2 x . . . Zpq x Z x . . . x Z (2.T) 

bi— times 

where C = {xiXjX^Xj 1 } is the commutator subgroup of F(X); Z p is the cyclic group 
(x\x p ); Z is the infinite cyclic group (x\ ); pi > 1; p q ^ 0; Pi\pi+i- The number of Z 
factors is called the first Betti number, b\, and the integers (j>i,P2, ■ ■ ■ ,Pq) are the torsion 
coefficients. 

If R is a square matrix such that det R ^ 0, then it determines a finite abelian group 
G of the order \G\ = ±detR. 

It is well known that not all groups can be fundamental groups of 3-manifolds and 
not all presentations of a given group can be realized in the form of a 2-skeleton of some 
3-complex. The obvious (but not the only one) source of the restrictions is the Poincare 
duality. Given a 3-complex C, the dual complex C is constructed by putting k- dimensional 
cells of C into correspondence to (3 — fc)-cells of C. A presentation of 7r(C) can be read 
off from the 2-skeleton, K 2 , of C. As we consider manifolds, tti(C) = 7r 1 (C). It implies 
a kind of duality between generators and relators in presentations of the fundamental 
groups of 3-manifolds. 

Let us fixe a spanning tree T of the dual 1-skeleton, K\. Relators in Eq. ( |2.1| ) are in 
correspondence with links of K\. This set of them is obviously excessive. A minimal set 
is given by links which is not in T e K\. In this case, relators of tti(C) correspond to 
generators of vri(C) and vice versa. In particular, their numbers equal each other. Such 
presentations are called balanced. 

3 Combinatorial meaning of the partition function 

Now we are in a position to interpret the partition function Z^(a,/<) introduced in 
Eq. ( P-.IQ ) in terms of combinatorial group theory. Every simplicial ball from deter- 
mines a set of ^3 + 1 generators of 7r 1 (C). Any 3-manyfold M 3 can be obtained out of 
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some ball by gluing its boundary triangles pairwise. This process induces corresponding 
identifications of edges and vertices. The glued triangles form a 2-dimensional spine^, 
K 2 C M 3 , having the Euler character x{^2) = 1- Every edge in K 2 gives a relator of the 
fundamental group, 7r 1 (M 3 ), determined by a cyclic order of oriented triangles attached 
to it. Thus we find 

Z«(«, / i)= Y, E e aN ^ a+ ^ (3.1) 

N 3 ,Ni {P N3+1 , Nl } 

where J2{p N +1 N } is the sum over all nonreduced presentations on A^ + 1 generators 
and Ni relators associated with the 2-spines of simplicial complexes from a given class 
([n 1 ] = {C: miC) fixed}): 

P n , m = (X n \R m ) = (oio 2 • • • an|rir 2 ...r m ) E [m] (3.2) 

where r^s are nonreduced words in af 1 . Every generator a { e X n appeares exactly 3 times 
in R m , simply because the triangle has 3 edges. 

Any spanning tree T of the 1-skeleton K\ of the 2-spine K 2 (T C K\ C K 2 ) fixes a 
minimal subset of the relators R'n 3 +i C -Rjvi- The deformation retract of a tree is a point: 
T \ *, and the fundamental group of any 3- manifold possesses a balanced presentation. 

It is convenient to weight every complex with the number of spanning trees T C K\ C 
K 2 . It can be done by introducing another system of free matter fields attached to vertices 
of triangulations. Then, analogously to Eq. (|1.10|) , we define the partition function 

Z( 2 )(«, (i) = Y det 'L det 'L e aNo ~^ = £ e"^ 3 E r ( H N s , a) (3.3) 

{C} iV 3 {H]Y 3 eM} 

where Z^j^ eM} S oes over a ^ balanced presentations from a given homotopy class: 
{Hn 3 = (Xn 3+ i\R' N3+1 )} = [7Ti]. All such presentations can be enumerated. Then, 
there exists an algorithm [IIJ allowing, in principle, for determining if a corresponding 
2-complex is a spine of some 3-manifold. 

T{H N „ a ) = Y, E e* = ^e*T No (^) (3.4) 

N {Pn 3 +i,n 3 +n ^H N3 } N 

where J2{p N3+1 N3+Nq ^h N3 } is the sum over all nonreduced presentations (Xn 3+ i\Rn 3 +n ) 
which can be deduced from a given balanced presentation H Ns = (X Na+ i\R N3+ i) following 
a pattern of a spanning tree in some simplicial complex. The last equality is a definition 
of T No (H N3 ). 

In this way we divide the counting problem into two steps: we estimate first the number 
of the balanced presentations in a class {Hn 3 } and then the number of simplicial complexes 

" By definition, a spine is a 2-skeleton in a cell decomposition with a single 3-cell. 
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giving any given balanced presentation from {Hn 3 }- The last number, Tn (Hn 3 ), is 
exponentially bounded because any configuration has the spanning tree structure and 
the bound is provided by the number of corresponding trees. To estimate the first sum, 
\Hn 3 G [tti]|, is a fundamentally difficult combinatorial problem. 

The balanced presentations enjoy the duality property. Therefore J2{h N3 } can be 
equivalently described as the sum over the balanced presentations in which the maximal 
length of relators does not exceed 3: £(rj) < 3, Vj. In this dual description, the cyllables 
ciid^ 1 and of never appear in the relators, because in that case there would exist a circularly 
closed edge, which is forbidden by the definition of the simplicial complex. Therefore we 
can consider reduced relators without loss of generality. 



4 Estimates for abelian presentations 

Let us formulate the following simple fact as 

Proposition 1. If the maximum length of relators equals 3 [£(rj) < 3, Vj], then there 
is at most 2 n balanced group presentations on n relators, (X n \R n ), having the same 
abelianization. 

Proof. As C{r,^ < 3, commutators, a6a _1 6 _1 , can not appear. Given 3 letters a, h and c, 
there are only 2 inequivalent relators having the abelianization abc: 

[abc] = [bca] = [cab] [acb] = [cba] = [bac] □ 

Let us consider the (non-trivial) case when vri(C) has the trivial abelianization: i?i(C, Z) = 
0. Balanced presentations of the trivial abelian group are in correspondence with the uni- 
modular integral matrices^ 

n 

P E = (x 1 . . .Xnl^RyXj) ReSL(n,Z) 
We are interested in the subset corresponding to triangulations: 

n 

M n = {RE SL{n, R) : £ R% < 3, Vi} (4.2) 

For all matrices R G M n , the allowed matrix elements are Rij = 0, ±1, and there 
are at most 3 non-zero elements in each row. We can independently permute rows and 
columns and multiply them by —1. It means the double coset 0(n, Z)\J\4 n /0(n, Z). The 
group 0(n, Z) = S n * Z 2 is the free product of the symmetric group S n and the group 
generated by the multiplication by the diagonal matrix diag(—l, +1, . . . , +1). 

w In principle, det R — ±1. We choose +1 for convenience. 



(4.1) 



11 



An R matrix can be further reduced. If the i'th row has only one non-zero entry, 
Rij 7^ 0, then we can remove it along with the j'th column by applying the AC2 and AC4 
moves. If there are two non-zero entries, Rij ^ and Rik ^ 0, we nullify one of them 
by using the AC2 move on columns. After that we can repeat the previous step in order 
to decrease the size of the matrix. Therefore, without loss of generality, we assume that 
(RR T )jj = 3 (j = 1, . . . , n) and {R T R)jj > 3 (j = 1, . . . , n). 

Proposition 2. Let a simplicial complex C give (maybe after the described reduction) a 
presentation matrix such that (RR T )jj = 3 Vj. If C is a 3-manifold, then J2i=i — 3, 
Vj. 

Proof. The reduced complex has only one 0-cell a , n 1-cells o\ and n 2-cells a\. Let us 
consider the regular neighborhood of its 1-skeleton K\ (cf. Ref. [[HI). The neigborhood of 
<r° is a 3-ball /3 3 . Its boundary is the 2-sphere df3 3 = S 2 . The neigborhoods of 1-cells a\ 
are solid tori. They intersect the sphere, df3 3 , in a disjoint collection of disks. The 2-cells 
a\ intersect d(3 3 in a number of simple arcs connecting the disks pairwise. The number of 
disks equals 2n, while that of arcs equals 3n. Any such configuration defines a graph in a 
natural way. If any edge E of a graph connects distinct vertices and the valency of each 
vertex V is at least 3, then Ij^E > 3#V. As in the case in question Ij^E = 3#V = 6n, 
we conclude that the valency of each vertex equals 3. □ 

Let us denote this restricted set of matrices as Ai^ = {R G SL(n, Z) \ (R T R)jj = 
(RR T )jj = 3 Vj}. If R has a block form, det R equals the product of determinants of the 
blocks. Let us consider the irreducible case. Any matrix R G M.^ can be transformed 
into the form R = I — B. Then B can be uniquely represented as the sum of two 0(n, Z) 
matrices B = P\ + Pi'- -Pi, 2 £ 0(n,Z). This form has the residual symmetry under 
conjugate permutations: R — > aRa^ 1 , a G 0(n, Z). Therefore we can take the conjugacy 
classes of the B matrices as representatives for the double coset. For example, we can 
specify the cycle structure of P{ 1 P2 and an action of P2 on it. Therefore, if the constraint 
det R — 1 is dropped, the total number of configurations grows factorially with n. The 
growth rate of the number of matrices giving det R = is factorial as well, because it is 
sufficient to equate any two rows or columns to vanish the determinant. 

A convenient pictorial representation can be given in terms of oriented graphs on n 
vertices, Q n . For it we interpret the B matrices as the adjacency matrices of the graphs. 
If Bij ^ 0, we draw an arrow from the j'th to the i'th vertex. There are exactly two links 
going away from every vertex and exactly two arriving at it (one corresponding to Pi and 
the other to P2). The determinant of R can be rewritten as the sum over all possible 
collections of oriented self-avoiding loops {Li . . . L p } in the graph Q n 

det(/-B)=exp(-j:— ) = 1+ E IT B a ( 4 -3) 

fc=l {l 1 ...l p } k=l (ji)eL k 

where Yl(ji)eL k g° es over links constituting a loop L^. 

If we interpret B as the incidence matrix (i.e., if R^ 7^ 0, the z'th edge is incident 
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to the j'th vertex), then corresponding graphs are always collections of disjoint circles. 
These give the cycle structure of P± P%. 

Obviously det R — 1, if the orgraph Q n has no closed contours. It does not necessary 
mean that Q n is a tree, because the loops in Eq. ( |4.3| ) are all oriented, i.e., one is allowed 
to travel only in directions shown by the arrows. If we had the stronger constraint, 
2j=i \Bij\ < 1 (i.e., exactly one arrow goes from every vertex), then the only possible 
values of det R would be 0, 1 and 2. In the case det R = 1, trees only would contribute. 
This simple model gives a solution to the analogous counting problem for 2-dimensional 
spherical triangulations (cf. Ref. @). 

Let us formulate the next intermediate result as 

Proposition 3. Let \JA^\ be the number of the equivalence classes of matrices {M e 
SL{n,Z) | M' = M", if M' = ^M"^ (u h2 e 0(n,Z))} obeying the additional con- 
straints (M T M) kk = (MM T ) kk = p, Vfe. Then, asn^oo and p finite, \M%\ < exp(An) 
with some finite constant A. 

Proof To prove that, we can use the following facts from the theory of group lattices 
(see, for example, Ref. 

(i) SL(n, Z) is a lattice in SL(n, R), i.e., SL(n, R)/SL(n, Z) is connected and has a finite 
volume with respect to any Haar measure on SL(n, R). 

(ii) A fundamental domain of SL(n, Z) in SL(n, R) can be conveniently described in terms 
of the Iwasawa decomposition: M = ujarj, where lj G SO(n); a = diag(ai,a 2 , ■ ■ ■ ,a n ) 
(with all a k > and Y[?=i a % — 1); V ^ N, where N is the group of upper-triangular 
unimodular matrices (i.e., if k > I, r\ k i = 0; on the diagonal r\ kk = 1; and, for k < I, r) k i 
is arbitrary). A Siegel set in SL(n,R) is a set of the form S tu = QA t N u defined by the 
inequalities: A t = {a : a, < ta i+ \ (1 < i < n)}, N u = {q : \r] k i\ < u (1 < k < I < n)} 
and Q coincides with the total SO(n) group. The fundamental domain is the Siegel set 
Stu with t = -7= and u=\. The Haar measure in this parametrization is 

dp, — duj JJ — JJ da.i Y[ drjij (4.4) 

i<j a 3 t=l i<j 

From the constraint 

«fc + Yl a ^ik = p ( 4 - 5 ) 

i<k 

we find that o? k < p for all k. Together with the unimodularity property of matrices, it 
gives 

n 

m i + Uvl)<P n - (4-6) 

3=2 i<j 

From the observation that 

(M T M)ij = afrjij + ]T a 2 k r} ki ri kj (i < j) (4.7) 

k<i 
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is a dilute matrix, it follows that rj nm is dilute as well: in any row of 77 there are only a 
finite number of non-zero elements. These equations mean that only a finite number of 
inequalities \rj k[ \ < ~ can be violated for every value of 

The action of SL(n, Z) divides the SL(n, R) group manifold into cells, each of them 
identical to the fundamental domain. The equation 

2 . , 

defines a hypersurface separating two different cells. 

Now, we can use a standard trick in analytic number theory: to estimate the number 
of distinct SL(n, Z) matrices obeying some given inequalities, we can simply count the 
number of cells inside a corresponding region in the SL(n, R) group manifold. 

As we are interested in counting the equivalence classes under arbitrary permutations 
of rows and columns, let us assume in addition that 

Oii > a i+1 . (4.9) 

For arbitrary k < n, 

f[a*= del , (M T M) 4J , (4.10) 

1=1 - J - 

and the permutations of rows rearrange the n-tuples changing values of a's, not simply 
permuting them. We can choose one in which the inequalities (^4.9| ) are satisfied. It is 
specific to our particular case. 

The unimodularity and the constraints (|4.5|) together with the inequalities ( |4.9|) imply 
that, in the n — > 00 limit, only a finite number of equations ( |4.8| ) can be simultaneously 
satisfied. If this number equals k, we can estimate the number of distinct configurations as 
(cn) k /k\ (with some constant c). As upper-triangular unimodular matrices form a group 
over integers, we can always obey the inequalities [77^] < 1/2 V(i,j). Analogously, for 
the left multiplications by SO(n, Z) matrices, we can choose any coordinates on the coset 
SO(n, R)/ SO(n, Z). Then, as follows from the theory of group lattices, the integrality 
constraints on matrix elements can have only one or no solution for any fixed tuple 
(«i, a 2 , ■ ■ ■ , Oi n ). 

The number of equivalence classes of dilute upper-triangular r)y matrices is expo- 
nentially bounded. The combinatorial reason for that is that they have, in a general 
position, more symmetry than densely filled matrices from the N group. If t] nrn = for 
k < n,m < I, then ar/a' 1 = 7] for all permutations a acting only on the rows and columns 
from k to / and identical on the others (a 6 l r < © S q © I r », q = I — k + 1). Let us take one 



14 



of such matrices and try to draw a consequence of this enhanced symmetry. We have 

a(ujar])a~ l = u'a'rf (4.11) 

As uj' = aucr' 1 G SO (n) and r( = ar/a' 1 G N, a' = acta' 1 has to represent an integer point 
in the a space (i.e., corresponding to a matrix from SL(n,Z)). Let £ min be a minimal 
distance between any integer points. If a, is permuted with etj and \ \og(ai/aj)\ « 
£ m i n , then, in a general position inside a lattice cell, there must be cui = aij, because 
otherwise there would be 2 different integer points in the fundamental domain. Thus, we 
conclude that, in a general position, the product ar] possesses the same symmetry under 
permutations as the dilute 77 matrix alone. Then, in the n — > 00 limit, the number of 
equivalence classes of such products has to be of the same order as for the 77 matrices. □ 

From the pure mathematician's point of view, the argumentation given above is not 
satisfactory. Hopefully, a mathematically minded reader will be able to complete the 
proof (or, maybe, find a loophole in it). Let us accept the claim of this proposition as a 
plausible hypothesis. Then we can conclude that 

Theorem 1 The number of 3-dimensional simplicial manifolds constructed of n tetrahe- 
dra and having the trivial homology group, H\(C,Z) = 0, grows at most exponentially 
with n: \C n \ < e Xn , for some finite constant A. 

Proof To finish the proof, we have to show that, for any matrix M G Ai^ from a 
given equivalence class, there are exponentially many simplicial complexes accociated 
with it. In terms of abelian presentations, the problem is reduced to estimating the 
number of inequivalent multiplications by upper triangular N 1 x (N 3 + 1) matrices. A 
multiplication by such a matrix from the left corresponds to a change of the basis of 
relators in a ~K\ presentation generated by AC2 and AC4 moves only. Let us remember 
that we are considering simplicial complexes. These changes of bases corresponds to 
formal 2- deformations of 2-spines of complexes. Therefore their number is restricted from 
above by the number of corresponding spanning trees in K\ C C. As the valency of the 
trees is equal to 4, the statement follows. □ 

Theorem 2 The number of 3-dimensional simplicial manifolds constructed of n tetra- 
hedra and having a fixed homology group Hi (C, Z) grows at most exponentially with n: 
\C n \ < e Xn , for some finite constant A. 

Proof As was pointed out in Section 2, an arbitrary presentation matrix can be trans- 
formed into a unique diagonal form by pre- and post-multiplications by SL(n, Z) matri- 
ces. Such multiplications correspond to changes of basis of relators as well as generators. 
Therefore, the statement is a consequence of the previous theorem. □ 

It is natural to assume that, having fixed the Betti number b\ only, one gets a factorial 
growth of the number of complexes as a function of a volume. Indeed, we can see that the 
number of inequivalent presentation matrices grows faster than an exponential in this case. 
On the other hand, it is a piece of mathematical folklore that there is always at least one 3- 
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manifold associated with a given abelian presentation^]. Although, such a manifold should 
not necessarily be simplicial, it is plausible assumption that the simpliciality condition is 
not restrictive enough to cut down the growth rate. 

5 Reduced model and connection with the 2-dimensional 
loop gas model 

As has been already discussed, in 3 dimensions there are several natural choices of classes 
of complexes to define the partition function over. And the question of universality classes 
is open. The simplest choice could be the set of 3-spheres with imbedded 2-collapsible 
spines. It looks very natural from the viewpoint of Section 3. We restrict acceptable 
presentations to the set of obviously trivial ones. Namely, to those that can be reduced to 
the empty one, (|), by using only AC2 and AC4 moves as follows from the little Wright's 
theorem. So, we arrive at the partition function for the reduced model (cf. Eq. (|3.3| ) and 

zS(«,/i)= E E E e^~^ = 

{K 2 cC\K 2 X*} {TCK2} 

= 2 Yl e aN °-» Ni (5.1) 

N3,No{P Ns+1 , N3+No ^(\)} 

where Y^(c^s 3 \ is the sum over all simplicial 3-spheres; Y^, 2 is the sum over 

1 ' ' {K 2 cC\K 2 /^*} 

2-collapsible spines, K 2 , embedded in C; J2{tck 2 ] is t ne sum over an spanning trees in 
K 2 . The abelianizations of all the presentations from the sum J2{p N3+1 . Ns+N ^(|)} have 
upper-triangular unimodular matrices. 

Let us remind a reader that, by definition, a spine is a 2-skeleton K 2 of a complex C 
in a decomposition with only one 3-cell. In the simplicial context, this 3-cell is a ball B3 
with a triangulated boundary. K 2 is obtained by some pair-wise identification (gluings) 
of all boundary triangles. If K 2 is 2-collapsible, one can use recursively only one gluing 
operation: the identification of 2 triangles sharing a common link on the boundary, in 
other words, a folding along the link. This can simulate any sequence of elementary 
2-collapses and expansions, with an elementary 2-ball being just a single triangle. 

In the dual language, one starts with an arbitrary spherical 3-valent fat graph and 
applies the move which can be represented as the flip of a link with the subsequent 

§We are indebted to V.Turaev for communicating this fact to us. 
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elimination of it: 



(5.2) 



Both of these moves are well known and have been used in Monte-Carlo simulations of 
triangulated surfaces [[/J. Having glued all triangles, one finishes with a collection of self- 
avoiding closed loops, the number of which equals Nq — 1. We weight every loop with the 
numerical factor e a . 

It is convenient, instead of erasing a links after a flip, to decorate it with the dashed 
propagator and keep its track in a process of subsequent foldings. To represent such 
configurations, we need to introduce the infinite set of vertices 



and so on 



(5.3) 



which are generated by flips. For example, 



gives 



gives 



and so on 



or 



or 



give 



or 



or 



or 



give 



and so on. 

In the end, we obtain planar diagrams with two types of propagators: solid ones 
produce closed loops while dashed form arbitrary clusters. The total number of vertices 
is equal to 2(N 3 + 1). 

Different sequenses of flips can lead to differently looking planar diagrams correspond- 
ing to the same 3-dimensional complex. The simplest example is given by the two possible 
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self-energy reductions: 



6 




It is easy to see that this overcounting exactly corresponds to taking a sum over all 
spanning trees {T C K2} in Eq. ( |5.1| ). More precisely, links chosen to fold along are not 
in T. 

Thus, the model can be represented as a non-local gas of self-avoiding closed loops 
with the partition function 

^(^) = E^^"^ (5.5) 
{V} 

where J2{v} is the sum over all the diagrams. The weight w-d is equal to the number of 
inequivalent starting configurations giving the same planar diagram T> (can be 0). The 
number of closed loops is simply equal to iVo — 1 and a can be identified with the inverse 
Newton constant. 

Let us take an arbitrary planar graph with dashed propagators. In order to obtain 
all possible starting configurations, we have to flip back dashed links in all possible ways. 
For any given dashed cluster, a result is a subset of all planar 3-valent graphs. Hence, the 
number of them is exponentially bounded as a function of the number of links. Let us 
give a few illustrative examples. First, flips have to be performed in linear clusters giving 
all possible planar 3-valent trees doubled in a "mirror" : 
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and so on. The number of initial configurations for a linear claster of n dashed links, C n , 
is given by Catalan's numbers generated by the function 



71+ 1 



oo oo V JLi ^ I 

n=0 n=0 V J ZX 




(5.6) 



After that dashed propagators attached to the ends of the linear cluster have to be 
reattached in all possible ways to links of these "doubled trees" . The number of config- 
urations can be estimated as follows. First, we have to calculate the generating function 
for the number of possible attachments of dashed propagators to a solid line: 



It can be obtained from the equation 



(5.7) 



q(u, n) 



+ 



1 — /! 1 — V 



— U/jL 



q(u,n) + 1 



(5.8) 



which gives 



q(v,n) 



(1-^(1-^) 



'\ — v — jj) 2 — {y + + u 2 jj? 



(5.9) 



Second, we have to know the generating function for the number of all possible at- 
tachments of dashed propagators to a /c-legged tree on the left (or on the right but not 
on both sides): 



fc+i 



t(x, v) := £n{::::(V) **!/»; t(x, 0) = C(x) 



(5.10) 



The equation for it is simply 



t(x,u)=xC(x) t -^- + l 
l — i/ 



(5.11) 
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which gives 



t(x, V) 



1 - V 



1 — v — xC(x) 



(5.12) 



Then, the quantity we need to calculate is 





T(x,u,n) = n{ 

n,m,k 



X V 



n,m,k 



k+1 




x k v n ^r 



x(l 



v) 



((1 - u) 2 - xC(xj) ((1 - /i) 2 - xC{x)) 



(5.13) 



Newly appearing clusters have to be expanded in all possible ways according to the 
same procedure. What is important is that their number is restricted by the coefficients 
of q(v,n), which is a factor in ( |5.13| ), and does not depend on a volume of the trees. 
Therefore, one simply iterates till all dashed links disappear. The combinatorial reason 
for the existence of an exponential bound for the number of obtained configurations is 
clear: they are all planar! 



We have just seen that the weights wt> in Eq. ( |5.5| ) are recursively calculable. However, 
it is difficult to really take account of the entropy of configurations and the model (|5.5|) 
seems to be analyticly unsolvable. Therefore, as was first proposed in Ref. ||, let us 
consider the localized version of it, namely, the dense phase of the self-avoiding loop gas 
matrix model |§ (in the sequel simply loop gas (LG) model): 



J LG 



d N2 Yd N2 Y[X u exp 



2 



Y exp tr Y 2 tr 2 log (l <g> 1 - \(Y ® 1 + 1 ® Y) 



(5.14) 



N ^2 n 

— tr Y 2 1 

2 2 

Here Y and X u are hermitian N x N matrices. We have attached the lower index to the 
gaussian X u variable to weight every closed loop with the factor n = e a . A = e _/1//2 . In 
the N —>■ oo limit, this model generates planar diagrams having only the simplest 3-valent 
vertices from the whole series ([5.3|) . For such diagramms, all the weights in Eq. ( |5.5|) are 
trivial: wv = 1- 
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Let us suppose that this truncated model bears some qualitative features of 3-dimensional 
simplicial gravity. It seems to be a reasonable assumption, because, as follows from the 
results of Ref. ||, the corresponding universality class is very large and includes all local 
perturbations of the model (|5.14[ ). On the other hand, multiple overlapping sequences of 
the flips can be regarded as a kind of renormalization group procedure, which could draw 
any given system to a stable fixed point in some imaginable space of all planar-graph 
models. It looks a bit more natural in dual terms. One takes an arbitrary spherical 
ball with a huge triangulated boundary and uses subsequently and randomly the folding 
operation (dual to the flip). Soon the boundary triangulation will be randomized. If 
this randomization has appropriate statistical properties, we could expect to find many 
models falling in the same universality class and, therefore, having the same continuum 
limit. 

The critical behavior of the loop gas matrix model is well known. For n < 2, it 
describes 2-d gravity interacting with c < 1 conformal matter 

c = 1 — 6- — — ; n = — 2cos7rg (5.15) 

9o 

while for n > 2 the corresponding matter is non-critical and the model trivializes. In this 
phase, the number of closed loops is proportional to the volume and the mean length of 
each remains finite. 

In the vicinity of a critical point /x c , the partition function behaves as 

Z LG ^(fi c -n) 2 ~^ (5.16) 
This formula defines the famous string susceptibility exponent 7 s t r . 
Two main quantities of interest are the mean lattice volume 

(N 3 ) LG = f- logZ LG « -(2 - 7 *r)-^ (5.17) 
and the mean lattice curvature 

(N )lg = f- log Z LG « (2 - 7str) n f ci ; n) - ^ log |/i c - u\ 
oa l^c{n) — /i oa 

*-^(N 3 ) LG + ^\og(N 3 ) LG (5.18) 
He oa 

where the coefficients can be calculated explicitly using the Gaudin and Rostov's exact 
solution [H 

9 11 d% tI 1 

— log/i c = -- — — ; —— = - — — (5.19) 

on 2 n + 2 on 2irgfismn(l — g ) 
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If 2-n< 1, then § 



(N )lg « + -^=10^)^ (5.20) 

4 7tv2 — n 

If n > 2, then 7 str is independent of n and (N ) LG is proportional to (N 3 ) LG . 

This behavior is strikingly similar to the results of numerical simulations. In Ref . Q , a 
phase transition in the model (|1.2|) with respect to the Newton coupling, a, was found. In 
the "elongated" phase, a > a c , the mean number of vertices, (N ), is strictly proportional 
to a volume, (N 3 ). In the "crumpled" phase, a < a c , (No) is a non-trivial function of (N 3 ). 
The analogy with the loop gas matrix model suggests that the most probable scaling is 
(N ) Rj Ci{N 3 ) + c 2 log(A^ 3 ) with c 2 singular at the critical point a c . Presumably, it is a 
type of this singularity that can be, in principle, calculated in continuum theory in order 
to compare predictions of both approaches. The obvious problem for such a hypothetical 
comparison is that a is a bare coupling. Therefore, in the lattice model, only critical 
points with respect to it could show some universal features. In any case, Eq. ( |5.20|) looks 
reasonable. For the dimensionful mean curvature, (R) = 2Tra(N ) LG + ca(N 3 ) LG , as a 
function of the voulume, V = cl 3 (N 3 )lg, we find 

(R)^^ 2 V + -^\ogV + ... (5.21) 
er v 2 — n 

As the scalar curvature has to undergo an additive renormalization, it is natural to find 
a linear volume term singular in the continuum limit: a — > 0, V finite. The logarithmic 
in volume term can have in this limit a finite coefficient, if simultaneously n — > 2 — 0(a 2 ). 
This is quite physical behavior. 

Of course, the loop gas matrix model cannot give us precise quantitative information 
about simplicial gravity. However, it is very plausible that it has qualitatively the same 
phase structure as models (|1.2|) and ( |3.3| ) and may be quite instructive from this point of 



view. 



6 Discussion 

1. Combinatorial group theory gives a natural mathematical framework and sets up a 
standard language for physical problems related to lattice models of 3-dimensional 
quantum gravity. All the formal group constructions with relators and generators 
have a natural geometrical realization in terms of 2-dimensional complexes (or fake 
surfaces, in a less formal parlance). And vise versa, geometrical constructions can 
be formalized in the group theory terms. It would be interesting to find physical 
models which could be formulated and solved entirely in terms of abelian presen- 
tations. It might be a mathematically adequate way to make physically meaninful 
approximations. 
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2. To determine a growth rate for different classes of complexes is the first necessary 
step in the process of constructing and investigating simplicial gravity models. It is 
an interesting field of research in its own right, altough more mathematical in spirit. 
Explicit forms of the asymptotic expansions would give rigorous and exact solutions 
of the correponding physical models. 



3. The status of the reduced model ( |5.1| ) is not yet absolutely clear. As it was formu- 
lated in Eq. ( |5.1| ), it looks as a genuine simplicial gravity model with some peculiar 
matter coupled to it. It is not clear if this matter has any local description in terms 
of fields in a manifold. The first natural question about the model to ask is if all 
simplicial spheres are really counted in Zj^(a,/i). Or, in other words, whether any 
simplicial 3-sphere allows for embedding at least one 2-collapsible spine in it. At first 
thought, Zeeman's Dunce Hat could serve as a source of counterexamples. However, 
triangles lying inside an initial ball can form such a fake surface. Even in the sim- 
plest case of 2 tetrahedra a Zeeman's Dunce Hat can easily be obtained. And, say, 



the complex constructed in Ref. [19| cannot be a counterexample because it is not 
simplicial. It seems to be an uneasy question in general. Say, the similar statement 
in case of simply-connected manifolds would automaticly prove the famous Poincare 
hypothesis. 



4. The model ( |5.1| ) can be successfully simulated numerically with help of the Monte 
Carlo technique. Results of the simulations as well as an account of the algorithm 
will be reported elsewhere PU[. Let us simply mention here that they support the 



hypothesis that both the reduced 3-dimensional model (|5.1|) and the 2-dimensional 
loop gas model belong to the same universality class in the continuum limit, with 
the Newton coupling playing a role analogous to the central charge in 2D gravity. 



It would be extremely interesting if the original model (|3.3j ) had a similar contin- 
uum limit. In principle, it would not contradict to what is known from numerical 
simulations of the pure 3D gravity model. If we remember that the analogous 2- 
dimensional matrix model discribes the c = — 2 matter coupled to gravity, which is 
topological, then it might seem less surprising that our 3-dimensional construction 
is equivalent to some 2-dimensional model. 
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